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Abstract 



Non topological solitons, Q-ba. 
As was shown by Cohen et al. 



Is can arise in many particle theories with U{\) global symmetries. 



a, 



if the corresponding scalar field couples to massless fermions, 
large Q-balls are unstable and evaporate, producing a fermion flux proportional to the Q-ball's 
surface. In this paper we analyse Q-ball instabilities as a function of Q-ball size ans fermion mass. 
\^ . In particular, we construct an exact quantum-mechanical description of the evaporating Q-ball. 

> ; 

00 This new construction provides an alternative method to compute Q-Ball's evaporation rates. We 

' shall also find the new expression for the upper bound on evaporation as a function of the produced 



fermion mass and study the effects of Q ball's size on particle production. 
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I. INTRODUCTION 



A scalar field theory with an unbroken continuous global symmetry admits a remarkable 
class of solutions, non-topological solitons or Q-Balls. These solutions are spherically sym- 
metric non-dissipative solutions to the classical field equations . In a certain way they 
can be viewed as a sort of Bose- Einstein condensate of "classical" scalars. The construction 
of these solutions uses the fact that they are absolute minima of the energy for a fixed value 
of the conserved f/(l)-charge Q. So in the sector of fixed charge the Q-Ball solution is the 
ground state and all its stability properties are due to charge conservation. An important 
amount of work has been done on Q-Ball dynamics and on their stability versus decay into 
scalars . Apart some existence theorems that depend on the type of symmetry and the 
potentials involved the stability of Q-Balls is due to the fact that their mass is smaller 
than the mass of a collection of scalars. 

In realistic theories the scalar field has a coupling with fermionic fields. The addition of 
this coupling modifies the criterions of Q-Ball stability since they can now produce fermions. 
This fact will have an important interest for cosmology since Q-Balls can play the role of 
dark matter 5,0]. Particle production from the Q-Ball will reduce its charge Q and at a 
certain point the Q-Ball will become unstable versus decay into scalars to finally disappear. 
This problem has been considered in 0, for the production of massless fermions by a 
large Q-ball. The method used was to construct the quantum field as a superposition, 
with operator valued coefficients, of the classical solutions. In most cases we can express 
the general solution as a superposition of partial waves. The next step will be to use the 
asymptotic behaviour of the fields, considering the far past and the far future where the 
fields can be identified to free ones. In most configurations the fields in the far past and in 
the far future are free fields and the relation, the S'-matrix, linking them together contains 
all the information needed to answer the question of particle production. The problem we 
have here is that the Q-ball is a time-dependent configuration, so we need to be careful 
when we identify the asymptotic states to free (static) ones. We must make sure that the 
identification is made before the interaction is turned on and after it is turned off. 

To avoid the problems linked to the time dependence of the Q-ball we shall propose 
an alternative method. This uses Heisenberg's picture of quantum mechanics. We shall 
construct the time independent state representing particle production, it is done by solving 
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the condition that no fluxes are coming from infinity (no particles are moving towards the Q- 
ball) . We can then build the Heisenberg field operator containing all the relevant information 
and time dependence. Particle creation is then computed by using the number operator, 
but any other operator valued quantity can be calculated. The use of this method needs 
no limit calculations on the Q-ball's size, so it can be used to study small Q-ball as well. 
The difficulty now lies in solving the production condition. This condition can be solved 
considering asymptotics of the fields far away from the Q-ball. 

The major difference between the standard construction and this one lies in the kinemat- 
ical conditions used. The standard S-matrix based method will solve matching conditions to 
compute all the reflection and transmission amplitudes of an incident wave. The approach 
used here is different in the sense that wc construct a state having no incoming wave, making 
the computation of scattered amplitudes useless. Using this construction implies that all 
the particles are created inside the Q-ball, since they all move away from the Q-ball. 

The other important fact to investigate is the production of massive fermions. Can a 
Q-Ball produce any type of fermions, and is the fermion mass relevant for the Q-Ball's fife 
time? To answer this question we can use both pictures, the problem is now that solutions to 
massive flcld equations have twice as many degrees of freedom. Even if the construction of 
the Heisenberg fleld operator is possible and not very difficult, the resolution of the particle 
production condition is a complex task to achieve, so we shall use the S'-matrix picture. 
The partial wave expansion of the solution can easily be obtained when working with one 
space dimension. In fact using this simple picture allowed us to obtain analytical results 
describing a fermion field being scattered by a Q-ball. Computing the evaporation rate as a 
function of the produced fermion mass will lead to a new definition for the absolute upper 
bound of evaporation rate. An other important question we can ask is the role of the Q 
ball's size on the different particle creation regimes. To answer this question we shall study 
both very big and small Q balls. 

The paper is organised as follows. We first give a review of the simplest 3 dimensional 
Q-Ball model to build up its basic properties and we then reduce it to 1 space dimension. 
We then consider interaction of Q-Balls with massless fermions and construct a solution 
where there are only fermions moving away from the Q-Ball, this construction will give 
us the evaporation rate of Q-Balls for production of massless fermions. We then consider 
production of massive fermions, using the incident wave outside the Q-Ball. This requires 
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the knowledge of the reflection and transmission amplitudes on the Q-Ball's surface. Finally 
we extend our results to 3 © 1 dimensions. 



II. A SIMPLE Q-BALL MODEL 

We review here the basic properties of a 3-dimensional Q-Ball using the simplest possible 
model. As we mentioned in the introduction, the Q-Ball is the ground state of a scalar theory 
containing a global symmetry. We can now build the simplest model in 3 © 1 dimensions 
having a Q-Ball solution: it is a S0{2) invariant theory of two real scalar fields (in fact it is 
the U{1) theory of one complex scalar field) 0. We start by writing down the Lagrangian 
and the equations of motion for the scalar field, to obtain the conserved charge and current. 
The Lagrangian of the scalar sector is given by : 

c = d,<j>*dy-um. (1) 

The U{1) symmetry is 



e*°c 



The conserved current is 



Jm = ^(0^9^0- (9^0^)0), (2) 



Q = / d^xjo. (3) 



and the conserved charge is 



We build a solution with the minimal energy : if f/(0) = is the absolute minimum of the 

aitential, = is the ground state and the U{1) symmetry is unbroken. It was shown in 
that new particles (Q-Balls) appear in the spectrum, if the potential is such that the 
minimum of is at some value 0o 7^ 0. 

Mm[2t//|0n = 2f/o/|0o|' < /i' = U"{0). (4) 
The charge and energy of a given field configuration are : 

Q = ±J(-dtr<P + c.c.)d^x, 

r • 1 

E = / i|0P + i|V0|2 + t/(0) d^X. 



The Q-Ball solution is a solution with minimum energy for a fixed charge, we thus introduce 
the following Lagrange multiplier 

e^ = E + uj[Q- ^J{(j)*dt^ + c.c.)d^x]. (6) 

Minimising this functional with the standard Q-Ball ansatz : 

= me'^\ (7) 

where 0(r) is a monotonically decreasing function of distance to the origin, and zero at 
infinity. Inserting the Q-Ball ansatz in the equations of motion gives in spherical coordinates 

If we interpret as a particle position and r as time this equation is similar to a Newtonian 
equation of motion for a particle of unit mass subject to viscous damping moving in the 
potential ^uj^(f)'^ — U. We are searching for a solution in which the particle starts at t = at 
some position 0(0), at rest, ^ = 0, and comes to rest at infinite time at = 0. Solving this 
problem is not difficult (see [l| for details). One of the solutions can be the localised step 
function. Although we can solve exactly the equation of motion, we will not do it in this 
work. 

This construction is the Q-Ball we where looking for, in the sense that it is the ground 
state of the theory with constant charge. We used only one field to describe it but it is in 
fact made of a collection of scalars. Q-Balls rotate with constant angular velocity in internal 
space and are spherically symmetric in position space. As the charge Q goes to infinity, lj 
approaches 



uo = V2t/o/02, (9) 

where Uq is the value of the potential at the minimum. In this limit, (p resembles a smoothed- 
out step function. The two regions {r < R and r > R) are separated by a transition zone 
of thickness This leads to the consideration of two approximations the thick and thin 
wall regime (see and jbl for details). The radius of the Q-Ball can easily be calculated 
using the definition of charge: 

Q = ^7tR'luo(I)1. (10) 



This calculation has been done 0(r) = if r > i?. All the properties of the Q-Ball are now 
known except the exact profile of the Q-Ball field. We shall now build up the Q-Ball solution 
to our problem. The energy is given by the integral (0) and using the previous properties 
and taking the limit V oo, where V is the volume of the Q-Ball, the energy becomes 

1 . 



E 



-UJ 



'V + UV. 



The charge becomes 



g = ^oi0iv. 



(11) 



(12) 



We wish to minimise the energy with fixed charge. Using eq. ()12j) to eliminate tu, we have 
in the limit of Q — oo, 

1 Q 



^ 2|0|2\/ 
As a function of V it has its minimum at 

Q 



V 



(13) 



(14) 



Here the energy is given now by 




(15) 



With some little modifications this construction can be adapted to any dimensions. The 
model we shall use in the next section is a 1 © 1 dimensional model. In one space dimension 
the Q-ball profile will be the standard step function localised in a space region of size /. 



III. PRODUCTION OF MASSLESS FERMIONS 



The solution to the problem of particle creation by a Q-ball can be solved using two 
different pictures. The first one is based on the S'-matrix formalism, using the idea that the 
field is free for t ±oo. This construction is done by finding the solution to the equations 
of motion for a fermion interacting with a Q-ball, in terms of a superposition of classical 
solutions. The quantisation is made by upgrading expansion coefficients to operators, this 
will give us the Heisenberg field operator. The S'-matrix will then be constructed by iden- 
tifying the fields in the far past and in the far future to fields having the exact positive 
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where the expansion was made using rotational eigenfunctions. Once the total solution is 
known it is simple to build the transformation from the far past to the far future. In the far 
past only the incoming wave will survive and in the far future only the outgoing ones. 

The construction we are going to use here is different, we shall in the first place solve 
the equations of motion and obtain the Heisenberg field operator representing a fermion 
interacting with a Q-ball. In one space dimension this solution will be expressed in the 
form, 



where the ipq^e, t, z) are a basis of the solution to the Dirac equation for fermions interacting 
with a Q-ball of charge Q. A{e) and B{e) are operators depending on energy, their anti- 
commutation relations are the standard ones if the ip solutions satisfy proper orthogonality 
conditions. The next step we shall use is consider the space asymptotics of this solution. 
Far away from the Q-ball {z = ±00 for one space dimension) the solution is the standard 
free field solution. This identification will give us a relation between the solution operators 
A{e), B{e) and the free asymptotic ones a{p), b{p). The only difficulty in this identification 
is that the quantisation of the solution was made using energy (due to the time dependence 
of interaction) while the asymptotical operators depend on momentum. The next step will 
be to define and solve the particle production condition, saying that no particles are moving 
towards the Q-ball. In terms of asymptotic operators it is 



The last step of the resolution will consist in using the total Heisenberg operator and 
the particle productive state to compute the fermionic fiux giving evaporation rate. The 
main idea used here is to construct a solution having no incoming wave, so we do not need 
to compute any refiected or transmitted coefficients. This way all the particles come from 
inside the Q-ball. As mentioned before this new construction gives a good alternative to the 
standard scattering kinematics. The other advantage of this picture is to allow us giving a 
consistent treatment of time continuity. 




aL(p)|\l/ >= 6l(p)|^ >= for p> 0, on the left 
ai{(p)|\E' >= bji{p)\'if >= for p < 0, on the right. 



(16) 
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In the next two subsections we shall build the solution and the relation to asymptotic 
operators, while in the two last subsections we shall solve the production condition and 
compute the evaporation rate. 



A. Solutions to the equations of motion 



Writing down the Lagrangian of a massless fermion having a Yukawa interaction with a 
scalar field gives in one spatial dimension, 

^ferm. = i^pa'^d^i/j + + h.c), 



(17) 



where the C superscript indicates the charge conjugated fermion. The equations of motion 
and their solutions are fully described in literature on the subject (|l[ O, O]). Instead of 
treating separately the fermion and the anti-fermion, we shall construct the exact global 
solution to this problem, this solution will be made of different parts first the solution inside 
the Q-Ball (for z G [— /, /]). Equations of motion for the two components of the \E' field are : 

{ido + id^)ipi - g(f)ip2 = 0, 



(18) 



and 



-iujQt 



in the zone from — / to +1 and zero everywhere else. Using the ansatz : 








V 







e 2 



/ 



A 
B 



(19) 



the equations of motion are reduced to the following 2x2 linear system 

^ ' 0. 



V 




k-e M 

M -(k + e) 

The determinant of the system gives k = ±-\/e^ — M"^ = ±k^. Solving for the two cases 
k = +ke and k = —k^, we obtain the solution inside the Q-Ball: 




A 



1 

M 




(20) 



where M = g(f)Q, g is the coupling constant and 0o the value of the scalar field. The second 
part is the solution when 0o = (outside the Q-Ball) it is, 

^ '11)1 \ I K ,1 \ 

. . (21) 



a 



L,R 
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where superscripts L, R indicate the left and right side of the Q-Ball. In order to solve 
Dirac's equation everywhere the solution needs to be continuous in space. Space continuity 
gives aX, z — —I : 

and at 2; = +Z, 

= Aa.e-^^^^-'^^ + Be+^('=^+^)'. 

These matching relations are used to express the solution only using the parameters coming 
from the inner part of the solution. This construction will allow us to build a state where 
there is no incoming fermion, all the fermions are now produced inside the Q-Ball. Putting 
together all these parts gives the full solution continuous in space and time : 




de 



\ 



z < -I 



J 



I 



2 /"g" 2 ' 



l>z>M 



where 



2: > +/ 



(22) 



OLe — 



M 



(23) 



A little work on the solution and on its orthogonality properties leads to a solution of the 
form : 



*Q = ^/ ^ee-* (^i,%{e)A{e)-rrQ{e)B{e)) e^^^Q(t), 



with 



il{t) 




(24) 



(25) 



and 



v(/2^(e,0)v 



Z < -I 



N± 



-ik^z 
-ikcZ 



ikcZ\ 



l> Z>+1 



V 



^ /±(e,-/)e-^ 

v(/2^(e,-ore-- 



z> +1 



(26) 



and the functions f^2 having the form 



me J) 



-ikel\ 



^ e*''(±a*e-^'=*^ + e*'). 



(27) 



Ar± = 47r (cosh [Im [fee] /](1 + \a, 
± cos[Re[/Ce]/]Re[ae]) , 



(28) 



Finally the time dependent matrix was introduced for simplicity, the N± are the normalisa- 
tion constants. Quantisation of solution ()24p was made using equal time anti-commutation 
relations for If the i/j"^ functions satisfy / dz{'ilj'^' {e')y'ijj'^{e) = 5a'aS{e' — e) we can show 
that, 

{A(e),At(e')} = JdzJ dz'{^^{z,e))\^^{z',e'))x{^^^ 



5{z'-z) 



The \E'q solution we obtained has now being upgraded to a Heisenberg field operator de- 
scribing fermions interacting with a Q-ball. 



B. Relation to asymptotic operators 



First we conjugate the second component of the above solution, in order to c omp are 
it with the standard free solution for a massless fermion in 1 © 1 dimensions (see for 
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details). We look at the asymptotic behaviour of the Q-Ball solution (j^ . On the left and 
right-hand side of the Q-Ball, it has to be the standard free solution since the interaction is 
zero outside the Q-Ball's volume. After elimination of integrals and standard manipulations 
and variable changes, we obtain at z — * — cxd : 

/ 



0{p) 
Oi-P) 



a{p) + 



e{-p) 



\ 0{p) 



bH-p) 



^ fi{e,l)B{e) 



+ 



ft{e,l)A{e) 



e=p 



e=p~ 



(29) 



6{p) is the heavy-side function, 6{p) = if p is negative and 6{p) ^ when p is positive. 

t / xt 



Multiplying eq. (|2g|) by and by 

oi-p) ) [ e{p) 



we obtain the two following equations 



'2n 



aiip) = [/+(e,/)A(e) + /r(e,/)i?(e)] 



"0 



0{p) + 



+ [f,^eJ)A^{e) + f,-{eJ)B^{e)] 



^bii-p) = [/+(e,/)A(e) + /r(e,/)i?(e)] 



e=p- 



+ [fneJ)A\e) + f,ieJ)B\e)] 



e=p+^ 



A-p), 

i-p) + 
e{p). 



(30) 



In these two equations the / subscript indicates we are on the left-hand side of the Q-Ball, 
the same manipulations on the right-hand side lead to : 
1 



/2tt 



anip) = [/+(6,-/)A(6) + /r(6,-/)i?(6)] 



e=p-^ 



+ [/+(e,-0At(e) + /2-(e,-/)i?t(e)] 



/ZTT 



0{p) + 
0{-p), 



(31) 



e=p 



^0 



+ [/2ne,-/)At(e) + /2-(e,-/)i?t(e)] 



e=p+- 



0{-p) + 

e{p). 



Checking the anti-commutation relations of operators ai^n and 6^, is a tedious task but 
using the different energy ranges and orthogonality properties of the f^2 it can be done. 
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These four equations will be the basis of the construction of particle productive state, since 
they give a relation between free operators (lower case) and solution operators (upper case). 
These relations will give the solution in terms of free operators, so the next task we need to 
achieve is to define and construct the particle productive state. 

C. Construction of particle productive state 

As mentioned before, the construction of this quantum state \1' will be done using the 
fact that there are no particles moving towards the Q-Ball. These are negative momentum 
particles on the left and positive momentum particles on the right. In terms of aifi, and 
hi^R operators : 



This construction will lead to the opposite sign of the fermionic current on the left and on 
the right hand side of Q-Ball using eqs. (jHOHSII)- We then obtain four equations. For positive 
p, we have : 



aL{j>)\^ >= bL{p)\'^ >= for p > 0, on the left 
a/j(p)|^ >= 6/j(p)|^ >= for p < 0, on the right. 



(32) 



(/+(e,/)A(e) + /r(e,/)i?(e)) 



^ >= 0, 



(/+(e,/))Mt(e) + (/r(e,/))^i?t(e) 



(33) 



^ >= 0, 



and for negative p 




• For e > +^ we only have the following two equations : 



(/+(e,/)A(e) + /r(e,/)fi(e))|M/>=0, 
((/2+(e, -/))M(e) + (/^-(e, -/))^5(e))|vl/ >= 0. 



(35) 



• For the negative range e < — ^ we have : 



((/+(e,/))Mt(e) + (/r(e,/))^i?t(e))|M/ >= 0, 
(/2+(6, -/)At(e) + f^{e, -l)B^{em >= 0- 



(36) 
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For the middle range e G [— ^7 + ^] we have : 

(/+(6,/)A(e) + /r(e,/)fi(6))|M/>=0, 



in 



-l)A^{e) + f,{e, -/)i?t(e))|vI/>=0. 



(37) 



(£0 

2 ■ 



(39) 



The range where |e| > +^ is easy to solve, since we expect the solution to be the vacuum 
and to lead to no evaporation. The determinant of the matrix : 

det[f ^'^''^^ ^'^''^^ ]] = 2{l-a% (38) 

[if2H^,-i)r (/2-(e,-o)v 

is always different from zero. The only solution for an evaporating state in this range is the 
trivial solution given by : 

A{e)\^ >= B{e)\^ >= for e > ^, 
At(e)|^ >= 5t(e)|^ >= for e< 

In fact these two equations are the same, because we can always use the transformation 
A{e) = A'{e)9{e) + B'\e)9{~e), all equations will have vacuum solutions. Here the anti- 
commutation relations are trivial to check because of the two different energy ranges. For 
the middle range e G [ 

where particle production occurs as first shown in 
range, only particles are created and changing the sign of Uq changes the particle type. We 
now need to check normalisation of these new operators describing the evaporating state 
and their anti-commutation relations. Defining the evaporation operators in all the energy 
ranges, we have 



+^] things are a little more complicated, this being the range 

2|. Taking a look at solution (j^UI) in this 



A^{t 



e < 



/8^(/+(e,/)A(e) + /r(6,05(e)) eG 



2 ' 2 J 



(40) 



A{e) e>+^ 



and 



bJe) 



2 



5t(e) e< 
Bit) e> 



6G[-f,+f] 



(41) 



2 



where the vStt factor is the normalisation 



The anti-commutation rela- 



tions of these operators are easy to check. They use the fact that \fi\^ = 4-- The particle 
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production state is now fully characterised by the simple relation : 

ae(e)|^>=6e(e)|^>=0. (42) 

This simple relation gives the ground state for a Q-Ball producing fermions. We could now 
compute lots of different properties of these fermions but we compute their number. The 
state defined by this relation has no incident fermion, it is exactly the state we wanted to 
build. The final step will now be to compute the fermionic flux and obtain the particle 
production rate. 



D. Particle production rate 



The particle production rate is given by the current operator j^{x) = iIj{x)'j'^iIj{x), which 
in our case is 



(43) 



that we shall apply on the evaporating state defined by the vacuum for Oe and operators. 
First we invert the systems pn]) and (PT|) to obtain : 



e > 



A{e) = alie) 
B{e) = blie) 

A{e) = ae(e) 
B(e) = 6e(e) 

A{e)- 
B(e) = 



87r2/+(e,0 



'Me) + blie)) 
Meie) - blie)) 



(44) 
(45) 
(46) 



VSn2f^ (e J) 

Now we can compute the first term of the current on the left hand side of the Q-Ball : Using 
anti-commutation relations and the separate range of integrals and the definition of Aie) 
and -B(e) in terms of evaporation operators ae(e), &e(e) we obtain : 



< o|^iVi|o > 



rfe(|/+(e,/)p < 0|a,(e)at(e)|0 > +|/r(e,OI' < 0|&e(e)6t(e)|0 >^ 



+ 



de 



/+(e,/) /r(e,0 



2r] 



2C 



< 0|6e(e)&I(e)|0 > . 



(47) 
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The other term of the current, proportional to 'ip24'2, is very similar but we need to be careful 
with the fact that ip2 is proportional to /^(e, I). Applying the same method we obtain : 



< 0|ae(e)aT(e)|0 > 



2/+(e,/) ' 2/r(e,/) — i-vw-evw,- - (4^ 
de (|(/2+(e,/))|' < 0|ae(6)at(e)|0 > +|(/-(e, 0)1' < 0\be{e)bl{e)\0 >) 

It is easy to check that |/2(e, Z)p = |/i(e, Z)p leading to the compensation of terms with 
infinite bounds. Finally the expression for the fermionic current on the left is : 



+00 



2 



JL 



2 

2 

2 

2 



de 



+ 



2/+(e,0 2/r(e,0 
cte sinhpi/i;^/] ^ 



(49) 



If the real part of fc^ equals zero (^ < M), we use the definitions 



K + e 
M ' 



Of 



and the current is then: 



jL 



de- 



sinh1-2v/M2 - 



In the limit MZ — 00 we can neglect the negative exponentials and the only factor we are 
left with is the | coming from the hyperbolic sine. The right hand side the current is the 
same except for the sign, so the total current is equal to one half. Using the continuity 
equation we can write : 

2jL- (51) 



djo , dji r^^dQ [ r. -i 



Jl-Jr 



The current does not depend on z so the value of the current is constant on both sides, so 
after integration of the current over e we obtain : 

dQ 1 



dt in 



(52) 



This expression gives the particle production rate as a function of Uq when uq is smaller 

n 

than M in the limit of big Ml. It is in fact ^2] an evaporation rate since it does not depend 
on the Q-Ball's size. The other importance of this result is that it gives a absolute upper 
bound on evaporation rate. The case when the imaginary part of is equal to zero is a bit 
more complicated to solve. In this case we need to explicitly compute the current integral, 
so this case will be studied numerically. The value obtained for the evaporation rate in ()52|) 
is the 1 © 1 equivalent of the results found in literature on the subject. 
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FIG. 1: Particle production rate for small values of Ad I and for a fixed value of = 0.5. 



1. Production rate in function of size 

We shall first consider the limit where / is small. In this case we write 
sinh'[2VM2 -e2/] = A{M^ - e^)P = A{Mlf{l - 



2(1-^^) = 2(1 -Mi)- (53) 



These two terms will simplify to give after integration over e : 



JL = I'M^^, (54) 

leading to the particle production rate: 

f^f^M^^... (55) 

This result enshures us the fact that when Ml = 0, the Q ball does not exist, the evaporation 
rate is zero. This behaviour is shown on figure ^ The next limit we shall study is the very 
large Q ball limit. To do so we take a look at the production rate for large values of the 
size parameter, Ml and observe that the production rate becomes constant for big values 
of the size parameter (see fig EI)- These considerations also stand for all the possible values 
of the frequency parameter the only difference is when ^ gets bigger the stability of the 
evaporation rate comes for bigger values of Ml. 
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FIG. 2: Particle production rate in function of Ml and for a fixed value of = 0.5 



2. Energy flux far away form the Q ball 



The next calculation we can do is the calculation of the energy flux far away from the 
Q-ball. In the case where we consider the observer very far from the Q-ball the only relevant 
coordinate is the distance to the Q-ball, we are in a one spacial dimension case. The energy 
flux a distant observer can measure is given after normalisation by M, 



dE 



.±2L ^ 

Mdtda J-^ ^M' 



sinh[2ifc(_^)^ 



2iki e -.1 

3 ^ M I 



_9 —2ik( 



(56) 



this expression is obtained by computing the energy flux through a sphere containing the 
Q-Ball. When the real part of equals zero the fraction becomes equal to one. The result 



in this range will be proportional to Uq 



3. Results of numerical integration 



We can now give the evaporation rate of a Q ball into massless fermions in function of 
its internal frequency. In the first figure we can observe a limit in the evaporation rate. The 
absolute upper bound can be computed using 

dN 



dt 



< 



"0 
2 



de = Uq, 



this absolute upper bound will be used to normalise the evaporation rate. 
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FIG. 3: Particle production rate "^Zflf as a function of ^rr in the limit of very big Ml parameter. 
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FIG. 4: Normalised particle production rate ^'^'^^ 
a very big Ml parameter. 



upper bound ^ function of ^ in the limit of 



IV. PRODUCTION OF MASSIVE FERMIONS 



The evaporation of a Q-ball into massive fermions is more complicated than the previous 
case. We can quite easily obtain the Heisenberg field operator but solving the evaporation 
condition is a difficult task, even with only one space dimension. So the method we are 
going to use is the same S*- matrix based method used in [3, O]- This picture will need as 
starting point the expression of the solution as a superposition of wave packets. It is done 
by expressing the motion equations in matrix form and then expanding the solutions over 
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the eigenfunctions. This gives the separation into left and right movers. Choosing which 
wave is the incident one, we can write the solution as 



(57) 
(58) 



where the u's and p's describe the solution away from the Q-ball and the L, R subscripts 
stand for the left- or right-hand side of the Q-ball. This solution has two incident waves 
associated with particles or anti-particlcs moving towards the Q-ball, giving two solutions. 
The reflected and transmitted waves are associated with particles moving away from the 
Q-ball. The same construction is done on the other side of the Q-ball, to give four solutions. 
Finally we obtain the total solution as a superposition of these four solutions with the 
expansion coefficients becoming operators. This canonical quantisation does not introduce 
any big problem and can be done in a straightforward way. The next step will be to 
consider that in the far past only the incoming wave survives, giving us a relation between 
the operators in the far past and in the far future (where only outgoing waves survives). 
The last step we shall do is compute the number operator. 

The only difficult task is the computation of reflection and transmission amplitudes ap- 
pearing in the solutions. Will shall provide two methods for calculating these amplitudes. 
One method will consist in calculating all the scalar products of the motion eigenvectors, 
while the other one will consist in the diagonalisation of the motion matrices. The results 
are fully consistent, and the two methods serve to illustrate a variety of physical insights. 
The main objective we shall reach is the computation of the new value of the upper bound 
on evaporation rate. This bound will now depend on the produced fermion mass and not 
only on the internal Q-ball frequency. 

The first three subsections describe the solutions in terms of eigenvectors of 4 x 4 matrices. 
In the three next subsections we compute the transmission and reflection amplitudes for 
massive waves on the Q-Ball's surface. In the last subsections we compute the evaporation 
rate and simplify the problem by diagonalisation of the motion matrices. 
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A. Preliminaries 



Using the same Lagrangian as for the massless case and adding a Dirac mass couphng 
for massive fermions, gives the fermionic Lagrangian : 



C = ■^i-^^d^.i) + g{i>^i)<^ + /i.e.) + Md{x)jP^. 



(59) 



(60) 



The equations of motion using two component ip-Held, in 1 © 1-dimensions are 

{idt + id,)7/ji - Me-^'^Va + MDip2 = 0, 
{idt - id,)ij2 + Me-'^V? + ^d^i = 0. 

These equations can be solved if we use fields with four degrees of freedom. Solving this 
system will give us the solution inside the Q-Ball, which is the first step. To solve these 
equations of motion we use the following ansatz 



^2 = gi{z)e<'-'^'>' + g2{z)e-'('+'^^'. 



(61) 



Due to the separation of time components this ansatz will lead to four equations. These 
equations can easily be modified to reduce the numbers of parameters : we divide all equa- 
tions by M 7^ 0. We shall now re-write these equations, taking 



fi{z) = Ae'P',f^{z) = Be'P%gi{z) = Ce'P'g*{z) = De'^'. 
After some re-arrangement of the equations we obtain : 



(62) 



(63) 



-92 + Mogi = pfi, 
e^gi - - Md/i = pgi, 

+ /i + Mof^ = pgl 

where e_ = e — ^ and e+ — e + This arrangement has the advantage that we can now 
write the ■0-field in terms of four component spinors as : 





















v 


{92 J 





(64) 
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The idea of using four-component spinors is that now the fermion field contains both energy 
components, just hke the solution used in the previous chapter. The other advantage is that 
this rearrangement leads to the standard four component spinor solution. The equations of 
motion become in matrix form, 
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— e_ 


Md 





-1 


\ 






(a\ 




-Md 


e_ 


-1 







B 




B 














= Mp 
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-Md 
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C 
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e+ 


/ 









(65) 



All the parameters are normalised by M. The fact that we have the M factor one the right 
hand side will allows us to simphfy the space components and replace I by Ml. All the 
parameters we are left with now are all dimensionless, we should read the matrix elements 
to be all divided by M and thus dimensionless. These satisfy 

tMt = M^, (66) 

with r = DiagonalMatrix[l, —1, 1, —1]. This symmetry will be used to perform the nor- 
malisation of eigenvectors. 

B. Solution inside the Q-Ball 

Using the dimensionless parameters and the eigenvectors of the motion matrix we can 
write the time independent solution inside the Q-Ball in the form : 

The v''s are the eigenvectors of the motion matrix while the p's are its eigenvalues. We shall 
not compute here the exact form of these eigenvectors, but we can prove that the eigenvalues 
are either purely imaginary or purely real. For reasons that will become clear later on, the 
first two terms of this solution have positive momentum while the two last have negative 
momentum. This arrangement does not modify the shape or any properties of the solution, 
but will greatly simplify the rest of the work. Inside the Q-ball the time dependent solution 
is : 

4 

^ = I] Cje'^'-^'K^fe'P^' + CJe-*(^+'^)* {v^°'^''e'P^'y , (68) 
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where the up superscript stands for the first two components of the eigenvectors, while the 
down one indicates we take the two last components. 



C. Solution without Q-Ball background. 
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Outside the Q-Ball the solution is given by the eigenvalues and eigenvectors of the fol- 
lowing matrix : 



(69) 



All its parameters are normalised by M , so they are all dimensionless. The eigenvalues are 
given by : 

Pi,3 = ±yjel -Ml^ ±pi, ^^^^ 

P2,3 = ±^4 -Ml = ±P2, 
Let us stress that all parameters are dimensionless since we have to read them as being 
divided by M, the Majorana mass couphng. The exact eigenvectors can be easily calculated 
but we do not need them. Here the pi,2 momentum can be complex or real. If we want some 
particle to propagate outside Q-ball we need both pi^2 to be real, it gives for e 



to solve this we must identify two cases. The first case is. 



e_ > ^ e > 



2 



e>Mn + '^> 
- 2 2 



e_ > Md 

the last inequality is verified if ^ > Md- The second case is. 



e-<0^.<f 
-e^>Mn^'^>^-Mo>e, 

once more the last inequahty is valid when ^ > Md- A similar calculation for |e+| > Md 
gives : 



e>MD 



UJQ 



and e < -Md 
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(71) 







FIG. 5: Sketch of the possible ranges for e 

The only way to avoid the gaps and have the two waves (both pi and P2) is for e to be in 
the range : 

ee[M,-^,^-M,]. (72) 

This range will mix both the particles and anti-particles and thus lead to the non-trivial 
Bogolubov transformation. We also have : 

f > Mo. (73) 

This range is the equivalent as the range defined for the massless case. Following the same 
construction as in the previous section, the static solution outside the Q-Ball can also be 
written in the form 

*o = EAV^"^^ (74) 
^o = Y.B,u,e^''^\ (75) 

this time the Bj coefficients are on the left-hand side of Q-Ball while the Aj are on the right 
hand side. The time dependent solution is once more given by, 

4 

These considerations allow us to find the particle production energy range where particles 
can propagate outside the Q-ball. The first part of our calculation is over, the next task is 
to compute the reflection and transmission amplitudes. It is in fact solving the matching 
equations in matrix form. The solution will give a relation from the far left to the far right 
of the Q-ball and the centre part (the Q-ball itself) will not appear directly. 
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D. Construction of scattering matrix 



We want to construct the matrix linking the solution at z = — oo to the solution at 
z = +00. We are searching for the matrix: 



(77) 



As a reminder the i?'s are on the left while the A's are the right hand side of the Q-Ball. This 
construction will not contain directly any inside parameters so the A and B can be considered 
as free parameters while the reflection and transmission amplitudes will be contained in the 
V matrix. To compute the matrix elements we need to solve the matching equations. 









B2 




A2 




= V 








^3 









E. Matching in space 

We first start by matching the solutions aX z = —I we have: 

Biup,e-'P'' + B2Up,e-'P'' + BsUp.e-'P'' + B^Up, 



g iPil — 



(78) 



We redefine the Bi and the Ci in the way 

Bi = BjQ ^' , Ci 



(79) 



Multiplying equation ()8()|1 by r, orthogonality property coming from the symmetry of the 
motion matrices : 

4 



i=i 



(80) 



Doing the same for all the 5's and writing down all the relations in matrix form we obtain : 



U 



(bA 




(cA 


B2 




C2 




= ss 




B3 













(81) 
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with, 



^ ll{ TVi u{ rt'2 u{ rf 3 uJtv4 ^ 
U2TV1 U2TV2 U2TV3 uT^TVi 

uJtvi uJtv2 uJtvs U^rVi 
y^u^TVi uJtv2 uJtvs uJtv4 J 



(82) 



Q-ipil 










-ip2l 

















(83) 



and 



U ^ 



( T 

'^pj'^pi 




V 














T 














We can then write for the expression we obtain at 2; = — Z 



B2 
B3 
^B,j 



U-^SE 



C2 
C3 

yC4y 



At z — +1 we have using the same definitions as for before : 



U 



A2 
As 
yA,j 



SE{-1) 



C2 
C3 



(84) 



(85) 



(86) 



Mixing up these two relations we obtain for the total transformation matrix V the following 
relation : 



^B^ 
B2 
B3 



= U-^S£E'~^S-^U 



A2 
As 



(87) 



25 



Using the definition 



s{i)s{-iy^ = E, 



(88) 



we can easily show that : 



[E,r] 



0. 



(89) 



As we shall find out later on the last form is a transformation that allows us to diagonalise 
the motion matrix inside the Q-ball. Using this matrix we shall construct all the refiection 
and diffusion coefficients for all the waves moving inside and outside of the Q-Ball. Before 
we continue we need to remember that = —pi and p4 = —p2 for both sets of p's (barred 
ones and no bar ones). We see here that the choice for normalisation of eigenvectors will 
just act on the U matrix that can be either the identity or the r matrix or any other choice 
we can make. Finally the diffusion matrix V we where looking for is given by : 



In fact nothing we shall do depends on these matrices but we have shown the full procedure 
for completeness. 

F. Construction of reflection and transmission amplitudes 

We shall first construct the refiection and transmission amplitudes from the left side to 
the right hand side of the Q-Ball. To do so we shall use the definition of the V matrix given 
by equation [77| The first two coefficients are linked to right-moving waves while the last 
two coefficients are linked to left-moving waves. This choice for arranging the waves was 
done for simplicity, we can show that any other arrangement leads to same results. In fact, 
we choose the simplest possible arrangement. Due to the shape of the u spinors and the Q 
matrix in front the first and the third coefficients of the free solution have the same energy 
while the second and the fourth coefficients correspond to another energy wave. We shall 
identify these two energy ranges to the type one particles (1) and type two particles (2). 



V = U-^SES~^U. 



(90) 
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Using equation [77| and separating the matrix into four two by two blocs we can write : 










^ Vn 


Vl2 


ri 




^V21 


V22 



\r2j 



(91) 



where the two arrows stand for the incoming waves, the first two coefficients will be replaced 
by one. Using the bloc separation of the matrix we find : 

(92) 
(93) 















(::) 


= V21 







leading to 




T 



(94) 
(95) 



= V2ivr/ 

n 

The TZ and T matrices will give the refiection and transmission amplitudes when they are 
applied on the incoming wave coefficients. These two matrices are two by two the first 
line corresponding to transmission or refiection of particles with two different incoming 
waves, while the second line gives the coefficients for anti-particles. We shall construct the 
transmission and refiection coefficients from the right to the left hand side of Q-Ball. Using 



the same method as before we 


have this time 
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(96) 



leading this time to 
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22 



(97) 
(98) 
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First case 



Second case 



FIG. 6: Sketch of both cases used to build the solution : we have each time two incident particles, 
two reflected and two transmitted. It is an effect of massive particles because we can not identify 
any more the particles with the anti-particles. 

Solving these two equations gives the reflection and transmission matrices for incoming 
particles from the left, they are : 

n = -Vi2Vn\ (99) 

r = V22 - V2iVri'Vi2. (100) 

Now that all the coefficients are known we can construct and quantise the solution. 
G. Construction of standard solution 

Using the transmission and reflection coefficients we can identify two different cases the 
first case is when incident particles are on the left hand side of the Q-Ball, while the other 
case stands for incident particles coming from the right hand side (see fig. E)). We shall treat 
separately the solution on the left and the solution on the right, the matching coefficients 
are those found in the previous section, they link the expansion coefficients on the right to 
those on the left. Writing down these two possibilities we have : 

= [Bie'^'^Up-, + B^e'^^^u,-, + ne-'^'^u.^-, + r^e-'^'^u.^.l (101) 

= [tie^^'^^Mp, + t2e*^"^"up-,], (102) 

for the first case, the two incident particles coming from the left hand side of the Q-Ball and 

= [tie-^^^^M-p, + t2e-*^"^"M_p-J, (103) 
= [ne'P^-Up, + he'^'^up, + A,e-'P^^u.p-, + A^e-'^'^'^u.^,] , (104) 

for the second case. In both of these definitions we have : 




(105) 



'n 



r2 




U 




(106) 



To clearly understand the construction, the 5's are the incident amplitudes from the left 
while the A's are the amplitudes from the right. The 1 and 2 subscript indicate the type 
of particle we are dealing with, we have two different exponentials in ^l{t). We then need 
to take the complex conjugate of the terms corresponding to the two last components of 
spinors. To continue building the solution we still need to separate each of these two cases 
in two, considering only one type of incident particle at the time. This construction leads 
to the four following pieces, that will be identified to the four degrees of freedom that our 
solution has : 



= [e^Pi^-Up-^ + me- 



ipix 



-Pi. 



for the one incident type one particle from the left and 



■P2\ 1 



-ipix 



U. 



■P2. 



+ [e*P2^Mp-2 + r22e" 



■tp2X 



u. 



-IP2X 



r2i 



[r2ie 

= hie'^^^'Up-^] + [t22e ""^^Mpaj, 

for an incident type two particle. The coefficients are given by : 

-f1 fi-f 

'^22 y \ 1 y \t22 J \ 

The two other pieces for particles incident from the right we have 



u. 



-P2ll 



*R = [riie^P^i^Mp, + e-*Pi"M3] + [ri2e-*^^^"«2] , 
for one incident type one particle from the right and 

*L = [i2ie-'P'^u-p-,] + [i22e-'P''u-p-,], 

= [f2ie'P'''Up-,] + [f22e-'P'''Up-, + e-'P^'^u-p-,], 

for an incident type two particle and finally the coefficients are given by 
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If we want to easily remember the coefficients there is an easy trick. For the coefficients 
without the tilde we read the subscript from the left to the right, r2i is the coefficient for an 
incident type two particle being reflected as a type one particle. The same lecture stands 
also for the tilde coefficients representing incident particles from the right. 



H. Quantisation and Bogolubov transformation 

Quantisation of solution is now easy, the total quantised solution will be a linear combina- 
tion of all four parts, with expansion coefficients becoming operators after the normalisation 
is made. The solution is given by : 

4 

^ = Y.{Aj, Bj)e'^'-^'^^u^Pe'P^' + {A*, B*)e-'^'+^^\u'^''e'P^'y , (115) 

in our case only the eigenvectors corresponding to the pi eigenvalue have up components, 
while only the eigenvectors corresponding to p2 have down components. 

= ei(^-^)t[eiP^^u-^ + me-'^'Vp-,] + e*(^+^)*[ri2e-^^^^xip-J , (116) 

incident particle is the wave containing Up^ , all this superposition must be represented using 
the same operator so we are sure to have only four degrees of freedom. Quantisation will be 
done using energy, for the incident wave we have : 

^ (117) 

-^e>MD + f, 

leading for this first wave to 



e^(^+^)*[ri2e-^^^2^(wl"p^")*]}6^(pi), (118) 

after the conjugation of the term proportional to e**^^"*"^-** we have, 

r ^ de{e^(^-^)*[e^^"i^ii^f + riie-^^i^ii!i^-j6j-^ + 

e-^(^+^ )* [r\^e^P> {u^^'Y] hp,}. (119) 

Applying the same method to all the terms we flnally obtain for the total solution having 
four degrees of freedom : 

poo 
J Mjj + -^ 
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J Mj^ — ^ 
poo 

f'OO 

J Mn — ^ 



The upper integration bound is ^ — Md so we are only left with the terms containing the 
a operators, 

J Mo — ^ 

(120) 

JMo—^ 

One interesting result can be found here is that like in the masslcss case if we change the 
sign of ^ we change the particle type since we change the operator type. For the moment 
the a coefficients are only expansion coefficients since we have not quantised the wave yet. 
The solution on the right hand side of the Q-Ball is given by : 



J Mn :r- 

(121) 

+ / t^e{[fi2e*^'i"Mp-JaLp-, + [f*2e-*P^X-2 +e'^'"«-p-2]«-P2} 



At i = +00 only the terms without any incident wave will survive we have then 



'^220(^2) + *22a(-P2) + ^210^(^2) + *i2a^(-P2) = aout{P2), (122) 

this is the Bogolubov transformation we where looking for. If we want the aout coefficient 
to be a operator we need to check that it satisfies the same anti-commutation relations as 
Qin- We have : 



{(aL)', aout} = ((r^i)'s- + (r*2)V4, + (4)'a_p. + (tt^Yalp,^) 

X ((^^12)42 + (^22)ap2 + (t2i)al 

P2 + (^22)cf-p2) 

= ((^2i)V2i + (r2*2)V22 + (4)'i2i + (4)'^22) {(afj', a,„}. 

This relation can also be obtained if we set that the incident current is equal to the outgoing 
one, or even with the normalisation of wave packets. At this stage it can be important to 
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use some normalised eigenvectors. As will shall show later on it is always the case if we 
diagonalise the matrix outside the Q-Ball. The number of created particles is now given by 

in < 0|aLt«o«t|0 >in= {-. |o , l^^^lo l-^^lo , 1 7 19 ) (^(e - e')- (123) 



We need to smooth out this result, to do so we shall use the same argument as [2] to finally 
obtain 
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de. (124) 



Since we are dealing with a Bogolubov transformation we have 

< 1, (125) 
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dN 

— <uo- 2Md, (126) 

In fact the best thing to do is to consider the solution in all space on the left and on the 
right instead of considering only one side. To do so we just need to consider an incident 
particle on the left and build the solution without tilde factors. The rest of the procedure 
is the same we have, 



^L+R = I " " deie^^^-'^^^ie-'^^'^ulM^ + e-^(^+^)*[e-^^^^^«^-")'^ + A^e'^> {ut"^;^r]a,,} , 

J M]j — Y 

(127) 

J M£, — ^ 

this time leading to 
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de, (128) 



dt 271 Jmd-^ 

after normalisation of operators. These results seem to be correct because when the fermions 
become massless there is identification of both types of produced particles so the total 
coefficient becomes equal to one as in the previous section and there is total reflection. We 
could stop our calculations here since we've got the expression of particle production state 
and all the matrix elements are known. But we can greatly simplify this problem so the 
expressions we obtain become simpler and more readable. 
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I. Direct construction of S Matrix 



Using the shape of the different matrices we deal with we think there is a simpler way 
to construct the diffusion matrix. In fact all the matrices of motion equations can be 
diagonalised using simple transformations that preserve the symmetry of the problem. If 
the matrices can be diagonalised the eigenvectors will have automatic orthogonality and 
normalisation properties. This diagonalisation will be done using simple transformations 
depending on six parameters. The symmetry to conserve is given in ()66|). 

Taking a look at the Mq matrix defined in eq. (jU^ we can diagonalise it using the Lorentz 
boost-transformation : 



(129) 



with, 



T 

rvi 



cosh(a;i) sinh(xl) 

sinh(xi) cosh(xi) 





^ 



cosh(a;2) — sinh(x2) 

- sinh(x2) cosh(x2) j 



T. 



(130) 



(131) 



The last equation ensures us the fact that tv'^t = v ^ and that the symmetry of the problem 



is conserved. Setting xi and X2 being solutions of 

cosh(2xi) 
cosh(2a;2) 

we find for the Mq matrix the following diagonal form 



e_ sinh(2xi ) 
e-l- sinh(22:2) 



where 



/C2 



/ci 

-A;i 

A;2 

-k2 

(M|,-ei)smh(a;i) 

{Ml-el)sinh(x2) 
Md 



(132) 



(133) 



(134) 
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All the parameters we find after this transformation are real, the fc's that we find represent 
the momentum of the particles in this new base. The same transformation applied on the 
Ml matrix defined by eq. (jH^ gives : 



sinh(x) 
- cosh(a;) 
k2 




- cosh(x" 
sinh(a;) 




\ 



(135) 



ki 

-ki 
sinh(x) cosh(x) 
^ cosh(x) sinh(a;) 

with X = Xi + X2- We can check that this new matrix has the same symmetry properties as 
Ml. This simple transformation allows us to eliminate the Dirac coupling of our equations 
but in presence of the Q-Ball it is replaced by a double Majorana coupling. Since this 
transformation is made every where, it will not change any of the properties. A solution of 
equations (|132p is easy to construct it is : 

cosh(2xi) 



sinh(2xi) 



coth(2xi) = 
2xi = argcoth( 



M, 



D 



(136) 



Once this transformation is made inside and outside of the Q-Ball the eigenvectors outside 
of the Q-Ball only contain one non zero component, the important fact here is that the 
Mq matrix is now self adjoint so its eigenvectors have the standard orthogonality properties 
without the r matrix. When we do the matching in space eq. ()8()|1 instead of multiplying by 
uJt we multiply by u\ and the S matrix is made of the components of M[ eigenvectors, so the 
only thing we did was to diagonalise the transformed matrix. To continue the Diagonalisation 
process we now transform M[ in the way : 



with, 



M'; = T{siS2fTM[{s^S2) 



cosh[?//2] sinh[y/2] ^ 

cosh[?//2] -sinh[y/2] 

-sinh[?//2] cosh[?//2] 

sinh[?//2] cosh[?//2] j 



(137) 



(138) 
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S2 



As before we have : 



cos[z/2] sin[z/2] ^ 

cos[z/2] -sm[z/2] 

-sm[z/2] cos[;2/2] 

sm[z/2] cos[z/2] ^ 



(139) 



(140) 



to preserve the symmetry of the problem. This set of transformations looks more complicated 
then the simple boost we used to start, it is the case for the parameters we shall need to 
use. But it is of great use for the final simplifications and results. Setting y and z to be 
solutions of : 

—2 cos{z) cosh(|/) sinh(x) 



sin(2;) = 
sinh(y) = 



ki - k2 
2 cosh(?y) cosh(a;) 

ki + k2 



(141) 
(142) 



we have 



with, 



^ tan(^) 
tanh(7/) 



—2 cosh(y) sinh(rc) 

ki - k2 
—2 cosh(a;) 
ki + k2 



A -M 

M -A 

B M 

-M -B 



A 1 //, , X / \ / 7 7 4 cosh(,T)^, , , , 
^ = o ^'1 - ^2 cosOz + /ci + A;2 - , , , ) cosh(t/) 
2 V ki + k2 



4 cos(2;) sinh^(,x) cosh^(t/) 
/ci - k2 



„ 1 // , , V /X /, , 4cosh(a;), , , , 
B^-[ {-ki + k2) cos(z) + {kr + k2- , , , ) cosh(?/) 

Z y Ki + K2 

4 cos(^) sinh^(a;) cosh^(i/) \ 



ki- k2 J 
cosh(?/) sinh(2a;) 



/ci + k2 



(143) 



(144) 



(145) 
(146) 
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Taking a look at this Ml matrix we see it has the same form as the Mq matrix so we shall 
diagonalise it using the same boost transformation. This time the transformation will not 
be a boost since some parameters can be complex. In fact before going any further we have 
to find the solution to equation (jl42|) that night be complex, ki + k2 is small so the fraction 
on the right hand side is always bigger than one. We shall have to set y = + rj, this little 
trick allows us to easily solve all these equations. Finally to finish the diagonalisation we 
transform using the Vi matrix : 

= TViTTV2tM[v2Vi 

= tvItM[v3, (147) 



with V3 = V2V1. This last transformation can also be done using a slightly different matrix 
the v[ matrix defined by : 



( 



cosh (a) sinh(a) 
sinh(a) cosh(a) 
cosh(6) — sinh(6) 



v 







— sinh(6) cosh(6) ^ 



T. 



(148) 



(149) 



We finally have for the M"' matrix the form 



Ml' 



^1 
0-^0 
6 

-6 y 



with 



where 



M^i = (^2 -M2)sinh(2a), 
Mi2 = (B^ - M2)sinh(26), 

M cosh [2a] = y4sinh(2a), 
M cosh [26] = Esinh(2a). 



(150) 



;i51) 



(152) 
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Using these transformation the diffusion matrix V can be expressed in terms of the diago- 
nahsation matrices in the way : 



This form will be in fact far more simple then all the other possible ones, so this is the reason 
why we decided to use it rather then the form with the scalar products of the eigenvectors. 
What we do is exactly the same since we work in a base where the matrix of motion equations 
is diagonal. If we keep the scalar products with the r matrix we see that the only to vectors 
having negative values are the negative moving ones, with a simple calculation we could link 
the T matrix to the helicity operator. This final expression we obtained will lead to simple 
results and is used to compute all the reflection and transmission amplitude the results are 
given in the next section. 

1. Small sized Q balls 

Using the results of previous section we where able to compute all the amplitudes for 
small sized Q balls. The method used was to replace the exponentials in the E matrix 
by : (1 — ipi.2,3,d)- These amplitudes still have complicated expressions but all of them 
exept ^22 are proportional to the size parameter Ml = I. In the limit where / goes to zero 
all the amplitudes fall to zero except ^22 going to one. The ^22 amplitude representing the 
probability of a fermion remaining a fermion. This probability is obviously one if the Q ball 
disappears. If the size parameter is small the amplitudes will become proportional as for 
massless particle production. This quadratic behaviour is shown on figure [7| 

V. RESULTS OF NUMERICAL INTEGRATION 

We first tested the stability of production rate in function of size to see if like in the 
previous case the particle production rate becomes constant and stable for big values of 
the size. If the production becomes constant above a certain size then we do not need to 
care about complex averaging processes. Figure |H1 shows the stability of evaporation rate 
for large Q balls. The little oscillations are due to numerical instabilities that vanish for 
very big values of size. We can now compute the evaporation rate for values of the Dirac 




(153) 
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FIG. 7: Particle production rate for small values of the size parameter. 
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FIG. 8: Particle production rate for small values of the size parameter. 



mass smaller then the Majorana mass (the coupling inside the Q ball). The next task we 
need to do is test the stability of our computations when the fermion mass parameter is 
bigger then the Majorana coupling inside the Q-Ball, it is the case when > 1. This case 
shows exactly the same behaviour of the other one except for the fact that it takes more 
computer time to obtain the plot of evaporation. The results are on figure ^1 A quick 
analysis of these results shows that there is a superior limit for all parameter sets, this limit 
does not depend on the mass parameter. Is seems to be normal since an infinite Q-Ball with 
an infinite internal frequency can produce any mass fermions. The last words we shall say 
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FIG. 9: Evaporation rate, for infinite (very big) Q-Balls in function of the frequency parameter 
for different values of the fermion mass parameter. 
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FIG. 10: Evaporation rate, for infinite (very big) Q-Balls in function of the frequency pa- 

rameter for different values of the fermion mass parameter bigger then one. 



about these results is that the "angle" in the curve correspond to the value of the frequency 
where the imaginary part of the impulsion inside the Q-Ball becomes zero, it is the point 
where particles start to propagate inside the Q-Ball. The normalisation of evaporation by 
its upper bound will lead to the same shape as the massless case but there will be a gap 
from zero to the value of the fermion mass. 

VI. ENERGY FLUX FAR AWAY FROM THE Q BALL 

The last step we need to achieve is compute the energy flux far away from the Q-Ball it 
is done by considering the flux through a sphere surrounding the Q ball. As before if the 
observer is far away from the Q ball the only important dimension is the distance to the Q 
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FIG. 11: Energy spectrum far away from the Q ball in function of ^ once for a very big Q ball 
and once for a small one. 



ball. We have, 

dE _ r+^-^ f |r2iP + |t2iP 

Mdtda ~ J-^+^a V|r2i|2+|r22P + 1^21^ + 1^22 

the transmission amplitudes disappear when the Q-Ball's size is very big. This integration 
can be done numerically and we can also introduce the Q ball's size to see its influence on 
the energy flux. The only difference is that a small Q ball will produce less energy for until 
the value of the frequency parameter becomes big. We could normalise these figures with 
the absolute upper bound. This normalisation does not introduce any new features since 
the normalised curve for very big Q ball would start with a constant part to then fall down. 
For the small one we will not find any constant part in the normalised curve. 
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VII. CONCLUSIONS 



As we have seen the couphng between the scalar field and fermionic field leads to par- 
ticle production from the Q-ball P|. To study this particle production we constructed the 
exact quantum-mechanical state describing the particle producing Q-ball. We used Heisen- 
berg's picture of quantum mechanics, the state describing the produced fermions is fully 
characterised by the fact that no fluxes are moving towards the Q-ball. This condition is 
solved considering the asymptotics of the flelds far away from the Q-ball. Using this state 
we constructed the Heisenberg fleld operator describing massless fermions produced by a 
Q-ball. This construction allowed us to prove that for large Q-balls in one space dimension 
the particle production does not depend on the Q-ball's size. While for small Q-balls the 
particle production rate is proportional to P. The extension of these results to three space 
dimensions is simple. For large Q-balls particle production is an evaporation, while for small 
Q-balls it depends on the size. The other result we need to point out is that we can consider 
a variety of kinematical constructions to compute the evaporation rates. The flrst one is the 
standard one where we compute the reflection and transmission amplitudes for an incoming 
wave. The second one we used is based on the fact that no particles are moving towards the 
Q-ball. We proved that these two pictures are equivalent. 

The fact that fermions acquire a Dirac mass does not introduce many changes. In 1 © 1 
dimensions the particle production rate does not depend on the Q-ball's size for sufficiently 
big ones. This result is not very surprising, since taking the limit m — * leads to the 
same results as the coupling with massless fermions. In this case the only difference is that 
evaporation occurs in a different range. The internal frequency uq, the energy of one single 
scalar forming the Q-ball, must be bigger than the produced fermion mass. This result is 
also quite intuitive, the scalars forming the Q-ball desintegrate into fermions so their energy 
must be bigger than the fermion mass. The second fact is that particle production occurs 
in the range mixing positive and negative frequency terms. In this range the Bogolubov 
transformations we build are non trivial. Using these two results we proved that evaporation 
can only take place in the range : [M/j — ^ — Md], with ^ > Mc- This result is in total 
accordance with the previous work done on the subject 0,0], and extends it a significant 
way. This new definition for the evaporation range will introduce a new upper bound for 
the evaporation rate. 
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When the Q ball's size becomes small there is no more evaporation since the production 
rate depends on the size. For small sized Q balls the particle production rate is proportional 
to The size will also slow down the energy flux a distant observer can measure. Taking 
the limit / oo does not need any complex averaging processes since the evaporation rate 
is constant and / independent for big values of the size. 

We also computed all the transmission and reflection coefficients for a massive fermion be- 
ing scattered by a large Q-ball. This construction allowed us to compute the exact profile of 
the evaporation rate. Using these profiles we proved that both constructions are equivalent. 
The last result we have proved is that evaporation rate is proportional to uq in one space 
dimension while it is proportional to uJq in three dimensions. In fact in three dimensions it 
is proportional to — M^)^. These reflection and transmission amplitudes will be used to 
study the behaviour of Q Balls in matter. 
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